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ABSTRACT 

A clear detection of excess of power, providing a substantial evidence for solar-like oscillations in 
the G5 subgiant /i Her, is presented. This star was observed over seven nights with the SARG echelle 
spectrograph operating with the 3.6-m Italian TNG Telescope, using an iodine absorption cell as a 
velocity reference. A clear excess of power centered at 1.2 mHz, with peak amplitudes of about 0.9ms~^ in 
the amplitude spectrum is present. Fitting the asymptotic relation to the power spectrum, a mode 
identification for the £ = 0, 1, 2, 3 modes in the frequency range 900 — 1600 /iHz is derived. The most 
likely value for the large separation turns out to be 56.5 /iHz, consistent with theoretical expectations. 
The mean amplitude per mode {I = 0,1) at peak power results to be 0.63 ms~^, almost three times 
larger than the solar one. 

Subject headings: stars: individual [pi Her) — stars: oscillations — techniques: radial velocities 



1. INTRODUCTION 

The search for solar-like oscillations in stars has ex- 
perienced a tremendous growth in recent years (See 
iBedding fc Kieldsenll2006l for a summary). Most of the 
results came from high-precision Doppler measurements 
using spectrograph such as CORALIE, HARPS, UCLES, 
UVES. In particular, recent measurements obtained with 
the SARG spectrograph have led to the determination 
of frequencies, mode ampl itudes, lifetime and granula- 
tion noise on Procyon A (jClaudi et al.l 120051 : iLeccia et al.l 
|2007() . a star for which the nature of oscillation spectruin 
is still debated (IChristensen-Dalsgaard fc KieldsenI 120041 : 
iMatthews et al.ll2004D . 

In this paper we report the detection of excess of power, 
providing evidence for oscillations, in the G5 subgiant 
pi Her (HR 6623, V = 3.417 ± 0.014, G5 IV), a slightly 
evolved solar- type star with mass 1.1 Mq, effect ive temper- 
ature Te = 5596±80 K, and log g = 3.93±0.10 (|Fuhrmannl 
I1998D . 

2. OBSERVATIONS AND DATA REDUCTION 

The observations were carried out over seven nights 
(2006 June 13-19) with the high resolution cross dispersed 
echelle spectrograph SARG (Gratton ct al. 2001) mounted 
on the Italian 3.6m telescope TNG at La Palma observ- 
ing station. Our spectra were obtained at i? = 144, 000 in 
the wavelength range between 462 and 792 nm, where the 
calibration iodine absorbing cell covers only the blue part 
of the spectrum (from 462 up to 620 nm). Exposure times 
were typically 60 s, with a dead-time of 55 s between expo- 
sures due to the readout time. The signal-to-noise ratio for 
most spectra was in the range from 200 to 400, depending 
on the seeing and extinction. In total, 1106 spectra were 
collected, with the following distribution over the seven 
nights: 27, 106, 183, 186, 180, 226, 198. The first two 
nights were affected by poor weather and technical prob- 



lems. 

The extraction of radial velocities from the e chelle spec- 
tra w as based on the method described by lEndl et aTl 
(|2000t) . The observed spectrum was fitted with a recon- 
structed one, by using a convolution between the over- 
sampled stellar template, the very high resolution iodine 
cell spectrum and the measured spectrograph instrumental 
profile. Essential to this process are the template spectra 
of /i Her taken with the iodine cell removed from the beam, 
and of the iodine cell itself superimposed on a rapidly ro- 
tating B-type star, which was the same for all the mea- 
surements. 

The resulting velocity measurements of /i Her are shown 
in the lower panel of Fig. [TJ They have been corrected to 
the solar system barycenter and no other correction, such 
as decorrelation or high-pass filtering has been applied. 
The rms scatter of these measurements is 2.53ms~^and 
the uncertainties for the velocity measurements were esti- 
mated from residuals in the radial velocity extraction pro- 
cedure and are shown in the upper panel of Figure. [l]where 
it can be noticed that most lie in the range 1-2 ms~^ , ex- 
cept for the second night where technical problems due to 
the guiding system occurred. 

3. TIME SERIES ANALYSIS 

The amplitude spectrum of the velocity time series 
was calculated as a weighted least-sa uares fit of sinusoids 
dFrandsen et al.' 1995; Arcntoft et"allll998l : IBedding etHI 
2 004 ; Kjeldsen et al. 2005) , with a weight being assigned 
to each point according to its uncertainty estimate ob- 
tained from the radial velocity measurement. We have 
then opti mized the weight following the approach dis- 
cussed in iButler et al.l |2004 This consists in (i) clean- 
ing all power at low frequencies (below 250 /iHz) from the 
time series, as well as all power from oscillations (800- 
1800 /iHz); and (ii) searching these residuals for points that 
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deviated from zero by more than it would be expected from 
their uncertainties. We found that 20 data points needed 
to be significantly down- weighted. 

Fig. 2 shows a pronounced excess of power in the power 
spectrum (PS) around 1.2 mHz which is the clear signature 
of p-modes oscillations for a G5 subgiant star. This feature 
is apparent in the power spectra of individual nights, and 
its frequency and amplitude are in agreement with theo- 
retical expectations, as we shall discuss. Moreover, it is 
clearly disentangled from the low frequency increase due 
to slow instrumental drifts (1// behavior in the amplitude 
spectrum) and the high frequency white noise. In partic- 
ular, we find that the mean noise level in the amplitude 
spectrum in the range 3-5 mHz is a = 0.1nis~^ which cor- 
respond to a mean noise level in the PS of 0.013 m^s^^ . 
Since this is based on 1106 measurements, we can deduce 
that the velocity accuracy on the corresponding timescales 
is 1.82ms^i . 

3.1. Large frequency spacing and oscillation frequencies 

The mode frequencies for low-degree, high radial order 
p-mode oscillations in Sun-like stars are reasonably well 
approximated by the asymptotic relation (Tassoul 1980): 

iy{n, I) = Aiy{n + ^ + e) - l{l + l)(5t'02/6 (1) 

where n and I are integers which define the radial order 
and angular degree of the mode, respectively; Az/ (the 
so-called large frequency separation) reflects the average 
stellar density, (5fo2 is sensitive to the sound speed gradi- 
ent near the core, and £ is a quantity of the order unity 
sensitive to the surface layers. Note that 6vo2 is the so- 
called small frequency separation between adjacent modes 
with / = and Z = 2. 

On attempting to find the peaks in our power spectrum 
matching the asymptotic relation, we were severely ham- 
pered by the single-site window function. As is well known, 
daily gaps in a time series produce aliases in the power 
spectrum at spacings ±11.57 /iHz and multiples, which are 
difficult to disentangle from the genuine peaks. Various 
methods for the search of regular scries of peaks have been 
discussed in the literature, such as autocorrelation, comb 
response and histograms of frequencies. We used the comb 
response function method where a comb response function 
CR(Ah') is calculat ed for all sensible values of (See 
iKieldsen et al.|[T995l for details) . In particular we used the 
following function: 





CR{Aiy) = Yl [PS{i 



2n-l 



'-Aiy)PSU^,,^±nAiy] 



(2) 

so that a peak in the CR at a particular value of Ai^ in- 
dicates the presence of a regular series of peaks in the PS, 
centered at z^max and having a spacing of Av 12. It differs 
from a correlation function as the product of individuals 
terms is considered rather than the sum. Moreo ver it gen- 
eralizes the CR defined in lKieldsen et al.l (|l995l ) including 
(and down-weighting) four outermost additional frequen- 
cies so that it better detects the presence of a regular series 
of peaks as we checked in simulated power spectra. 

To reduce the uncertainties due to the noise, only fre- 
quencies with amplitude greater than 4 ct in the amplitude 
spectrum, in the frequency range 800 — 1800 /zHz, have 



been used to compute the CR . We then determined the lo- 
cal maxima of the response function CR^Av) in the range 
40 < Av < 80 /iHz and the resulting cumulative comb 
response function, obtained by summing the contributions 
of all the response functions, had the most prominent peak 
at 57/iHz as shown in Fig. 3. 

The large separation Av scales approximately with 
the s quare root of the mean density of the star (|Coxl 
Il981h and by e xtrapolating from th e solar case {Av = 
134.8 /iHz; see iBedding et all 12001 one obtains Av = 
(M/Mo)i/2(i?/i?Q)-3/2 ^j,^ = 59 ± 4/iHz where the er- 
ror derives from the uncertainties on the radius quoted in 
iFuhrmaiml ()1998[ ). Since we searched over the entire range 
40-80 /iHz, the agreement between observation and theory 
is encouraging. 

We tried to extract the oscillation frequencies directly 
from the PS by using the modified extraction method 
described in IKieldsen et al.l (|1995D and recently used in 



iLeccia et al.l ( 2007( ): for each frequency ^-max extracted in 



the step before, a frequency region of 2 /iHz width centered 
on fniaxi'T^ 57 /iHz, with n = 1, ... was selected. When this 
region contained a peak with an amplitude greater or equal 
than 3 fj in the amplitude spectrum we identified this peak 
as the second component to be cleaned, and recomputed 
the PS for each n. The advantage of using this approach 
is that the extracted frequencies are determined except for 
the same day-night alias shift for each group of frequencies 
extracted during the procedure. In order to identify the 
individual frequencies of the modes we then constructed 
echelle diagrams corresponding to values of the large sep- 
aration around 57 /iHz, so that we could easily identify the 
ridges of the £ = and then those of £ = 1, 2, 3. 

It is important to see if there are additional frequen- 
cies, not following the comb-like structure, which where 
not found with the modified extraction. We thus analyzed 
the spectrum of the residuals obtained with the modi- 
fied extraction, as shown in the upper panel of Fig. 4, 
and we found 5 additional frequencies with amplitude 
> 3(7 namely, vi — 1133.6 /iHz; V2 ~ 1580.5 /iHz; 
V3 = 1387.1 /iHz; V4 = 1172.9 /iHz and V5 = 1045.6 /iHz. 
These frequencies can be of stellar origin but can also 
be spurious peaks due to the noise or due to the non- 
linearity of the extraction procedure. On the other hand 
it is possible to note that that V4 can be identified with 
an / = 3,n = 18 mode and v^ — 11.57 = 1034.03 can be 
identified with a Z = 2, n = 16 mode. The left panel of Ta- 
ble 1 shows the identified frequencies according with the 
modified extraction. It should be noticed that although we 
define a peak to be significant if its amplitude is greater 
or equal than 3 cr in the amplitude spectrum, in Table 1 
we have included one frequency with amplitude 2.9 cr as 
this frequency appears also with the standard extraction 
method (see below) and with an higher signal to noise ra- 
tio. 

The large and small frequency separation and the con- 
stant £ have thus been determined by means of a least 
square best fit with the asymptotic relation ([1]) for all 
the identified modes by means of the modified extrac- 
tion, which provided us with: Aiy = 56.46 ± 0.11 /iHz, 
Svq = 5.22 ± 0.81 /iHz and e = 1.46 ± 0.05. The expected 
position of the modes according to this solution is depicted 
with dashed lines in the left panel of Figure 5 while Fig- 
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ure 6 shows a close-up of the power spectrum with the 
identified frequencies according with the modified extrac- 
tion. 

As deviations from Eq. (1) coul d be expected for (. ^ Q 
in the case of avoided crossings (jAizenman et al.|[l977D . 
it is important to check that the standard iterative sine- 
wave fitting leads to a consistent identification. We thus 
extracted iteratively all the frequencies with amplitude 
greater than 3 a and tried to construct an echelle diagram 
by searching in the range 50 — 60 /iHz. It turned out that 
only with Az/ « 57 it was possible to clearly see the ridge 
corresponding to = 0. It was then straightforward to 
identify the £ = 1,2,3 mode, except for two frequencies, 
namely vi = 1112.0 /xHz, V2 = 1596.6 /zHz which could 
not be identified. The PS of the residuals in shown in the 
lower panel of Fig. 4 while the resulting echelle diagram is 
depicted in the right panel of Fig. 5. 

A synoptic view of the extracted frequencies with the 
two methods is then presented in Table 1 where the fre- 
quencies extracted according with the modified extraction 
(left) are compared with those obtained with the standard 
extraction (right). It is reassuring to notice that 11 fre- 
quencies have the same identification in both approaches. 
In particular the identification of the £ = ridge appears 
to be quite robust, and indeed one can see in Fig. 5 that 
the echelle diagrams obtained with the modified extrac- 
tion (left) and the standard extraction (right) are essen- 
tially the same; this result suggests that we have correctly 
detected the large separation in this star. 

In Table 1 the error on frequencies is reported in paren- 
theses. In particular we have used the analytical estimates 



CT(a) = Y <^{'") - 0-1 m/s, 

.(^) = ^ ^ - ^ .Hz (3) 
\ N ttT a a^/s 

()Montgomerv and O'DonoghudlTOQll ) where N = 1106 
is the number of data points, T = 6.25 is the total du- 
ration of the run in days and cr(u) = (ay) is the mean 
error of the data. Note that a (a) represents the noise of 
the amplitude spectrum after pre-withening and a (a) /a is 
the inverse of the S/N ratio. However the errors given 
by these formulae represent a lower li mit to the real error 
(|Montgomerv and O'Donoghud Il99ll ). In particular the 
actual errors can be several times larger up to a factor 
« 10 (note that the formal resolution is 1.8 /iHz). 

The large and small frequency separation and the con- 
stant e have thus been determined by means of a least 
square best fit with the asymptotic relation ^ for all 



the identified modes, which provided us with: Aj/ = 
56.50±0.07 /xHz, Si^o = 5.03±0.94 /LtHz ande = 1.44±0.03. 
The expected position of the modes according to this so- 
lution is depicted with dashed lines in Figure 4. Figure 5 
shows a close-up of the power spectrum with the identified 
frequencies according with the modified extraction. 

3.2. Oscillation amplitudes 

The strongest peaks in the amplitude spectrum of /j Her 
(square root of power) reach about 0.9 ms^^ . However, 
these are likely to have been strengthened significantly by 
con structive interference with noise peaks. As stressed 
by iKieldsen fc BeddinS (jl995D the effects of the noise 
must be taken into account when estimating the ampli- 
tude of the underlying signal. In order to estimate the 
average amplitude of the modes w e used the result of 
IKieldsen fc BeddinsI (|1995D (see also IKieldsen et al.|[2005l 
), where the effect of the noise is estimated. By con- 
sidering the five strongest peaks we found a mean am- 
plitude of the p-mode peak power, for the modes with 
I — 0, 1 in the frequency interval 0.8 — 1.6 mHz, of 
(63.14 ± 3.0)cms~^. The observed amplitudes are thus 
in good a greement with the t h eoreti cal amplitudes ob- 
tained by IKieldsen fc BeddinS (|l995f) for which Vosc = 
{L/Lq)/{M/Mq) (23.4 ± 1.4) cm s'^ = (69^ 50 ± 
4.1) cm s~^ (see also the discussion bv lHoudek1l2006( l. We 
also note that the frequency of excess of power (1.2 mHz) is 
consistent with the value expected from scaling the acous- 
tic cutoff from the solar case. 

4. CONCLUSION 

Our observations of Her show an evident excess of 
power in the PS, clearly disentangled from the low fre- 
quency increase, and with a position and amplitude that 
are in agreement with expectations. Although hampered 
by the single-site window, the comb analysis and the 
echelle diagram show clear evidence for regularity in the 
peaks at the spacing expected from asymptotic theory. 
Moreover, our results provide valuable confirmation that 
oscillations in solar-like stars really do have the ampli- 
tudes that scales as L/M by extrapolating from the Sun. 
We hope that in the near future a multi-site observing 
campaign will allow us to explore further the oscillation 
spectrum of fi Her. 

We would like to thank Tim Bedding and Antonio 
Frasca for important suggestions and encouragements. 
This work was supported financially by the INAF grant 
"PRIN - 2006" . 
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Table 1 

Oscillation frequencies in /i Her for the modified extraction (left) and standard extraction (right) 



(mHz) 


S/N 


Mode ID 




(mHz) 




S/N 


Mode ID 


985.5(0.2) 


5.3 


I = 0, 


n = 


16 


985.6(0.2) 




4.9 


1 = 0, n = 


16 


1100.4(0.3) 


3.8 


I = 0, 


n = 


18 












1211.2(0.1) 


9.1 


/ = 0, 


n = 


20 


1211.2(0.1) 




9.1 


1 = 0,71 = 


20 


1268.7(0.3) 


3.1 


I = 0, 


n = 


21 










1323.2(0.3) 


2.9 


I = 0, 


n = 


22 


1323.4(0.3) 




3.5 


l = 0,n = 


22 


1437.6(0.3) 


3.5 


I = 0, 


n = 


24 


1437.8(0.2) 


=1426.1+11.6 


5.3 


1 = 0, n = 


24 












1608.2(0.3) 


=1596.6+11.6 


3.3 


1 = 0, n = 


27 


1125.1(0.2) 


4.1 


I = 1, 


n = 


18 


1124.8(0.1) 


=1136.4-11.6 


6.9 


I = 1, n = 


18 


1180.0(0.2) 


4.2 


I = 1, 


n = 


19 


1180.0(0.2) 


=1192.4-11.6 


5.0 


I = 1, n = 


19 


1238.9(0.3) 


3.8 


I = 1, 


n = 


20 










1296.4(0.3) 


3.0 


I = 1, 


n = 


21 












1351.9(0.2) 


4.4 


I = 1, 


n = 


22 


1352.1(0.2) 




4.7 


I = 1, n = 


22 


1034.0(0.2) 


3.2 


I = 2, 


n = 


16 


1034.3(0.3) 


=1045.9-11.6 


3.3 


l = 2,n = 


16 












1205.6(0.3) 




3.2 


l = 2,n = 


19 












1375.4(0.3): 


= 1387.0-11.6 


3.4 


l = 2,n = 


22 


947.6(0.2) 


4.0 


1 = 3, 


n = 


14 


947.5(0.3) 




3.9 


I = 3, n = 


14 


1061.0(0.2) 


4.0 


1 = 3, 


n = 


16 


1061.2(0.3) 




3.9 


I = 3, n = 


16 


1173.0(0.3) 


3.2 


1 = 3, 


n = 


18 


1172.9(0.3) 




3.1 


I = 3, n = 


18 












1285.1(0.2) 




5.1 


I = 3, n = 


20 












1569.2(0.3) 




3.5 


I = 3, n = 


25 
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Fig. 1. — Velocity measurements of Her obtained with SARG (lower panel) and the corresponding uncertainties (upper panel). 
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Fig. 2. — Power spectrum of the SARG velocity measurement of // Her with the spectral window (in the inset). 
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Fig. 3. — CR function of the power spectrum of the n Her data. 
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Fig. 4. — The power spectrum of the residuals obtained with the modified extraction (upper panel) and with the standard extraction (lower 
panel) . The dot-dashed and the dashed lines indicate 4 a and 3 a respectively. 
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Fig. 5. — EchcUc diagram for the extracted frequencies for the modified extraction (right) and standard extraction (left). Squares, crosses, 
rhomboes, triangles are for £ = 0, 1, 2, 3 respectively. The dashed vertical lines represent the position of the € = 0, 1, 2, 3 modes according to 
our fit of the asymptotic relation. 
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Fig. 6. — Close-up of the spectrum of /^ Her with the identified frequencies. Dotted lines are the day-night aUases. 



